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BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to a method for modeling semiconductor devices 
and more particularly to a method for modeling semiconductordevices, such as field 
effect transistors (FET) and high electron mobility transistors (HEMT) which utilizes 
a semi-physical device model along with a data fitting model to form a relatively 
accurate physical device model for the use in large-signal microwave circuit computer 
aided design (CAD) tools. 

2. Description of the Prior Art 

HEMT technology provides unparalleled, high-performance characteristics at 
high frequencies (microwave to millimeter wave). As such, HEMTs are used in various 
RF applications. In order to accurately forecast the performance of such devices it is 
necessary accurately model the effect of the components physical structure on its high 
frequency, nonlinear and large signal characteristics. Thus, it is necessary to know how 
physical changes to the device will effect device performance in order to determine 
what process changes may be acceptable to improve RF yield product and which may 
be unacceptable which decrease RF product yield. 

Large signal microwave circuit computer aided design (CAD) is used to predict 
the performance of a modeled circuit at high microwave input power levels. The bulk 
of large signal microwave circuit CAD is based upon a harmonic balance method, 
which is a frequency domain analysis. This method demands that elements within the 
circuit satisfy sel f-consistent conditions, such as charge conservation for all steps within 
the method's algorithm. 

Currently, there are two predominant methods for modeling large signal 
characteristics of HEMT devices: table-based device modeling; and data-fitting non- 
linear device modeling. 

Table-based device modeling utilizes measured, bias-dependent small signal 
characteristics for a device along with software-based algorithms to represent a full 
large signal model for the device. Because the measured data usually spans almost the 
entire range of the possible device biases, a large table of data must be created. As 
such, this method is known as a table based or database model. Such table-based 
models are disclosed in: "A Simplified Broad-Band Large Signal Non-Quasi-Static 
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Table-Based FET Model," by M. Fernandez-Barciela, et al., IEEE Trans., Microwave 
Theory Tech , Vol. 48, No. 3, March 2000. The algorithms that make up this method 
are designed to produce modeling expressions that satisfy the charge conservative 
constraint required in a harmonic balance simulation. Thus, table-based method are 
known to perform well in large signal microwave circuit CAD. 

Although the table-based method has proven to be accurate and versatile for large 
signal microwave circuit CAD, it has several disadvantages. For one, the method 
requires the use of special, software-based algorithms that are only available through 
some commercial software packages. Some of these packages include the microwave 
design system (MDS), advanced-design system (ADS) and ICCAP by Agilent 
Technologies. Secondly, this method is not based on the device physics and thus does 
not describe the device's performance in terms of any of its physical characteristics. 
Thus, the method is not able to predict device performance based upon physical data. 
Microwave circuit designers and device technologists have developed empirical data- 
fitting analytical expressions that can both fit measured DC and small signal 
characteristics of semiconductor devices and fit measured DC and small signal 
characteristics of semiconductor devices which satisfy the charge conservation law. 
Thus, the so-called non-linear data-fitting device models perform robustly in harmonic 
balance based large signal simulations and are thus ideal for use for large signal 
microwave circuit CAD. However, because these models are completely empirical, 
they are not able to relate real, physical characteristics with the predicted device 
performance. Examples of data-fitting non-linear models for HEMT devices are 
disclosed in: "A New Empirical Non-Linear Model for HEMT and MESFET Devices," 
by I. Angelov, IEEE Trans. Microwave Tech. . Vol. 40, No. 12, December 1992. Thus, 
there is a need to provide a large signal microwave circuit CAD model that is able to 
relate physical device characteristics to its predicted performance. 

SUMMARY OF THE INVENTION 

Briefly, the present invention relates to a hybrid modeling method that forms 
models from a semi-physical device model along with an accurate data-fitting model 
in order to implement a relatively accurate physical device model as a model that is 
usable in large signal microwave circuit computer-aided design (CAD) tool. The semi- 
physical device model enables accurate representation of known physical device 
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characteristics and measured bias-dependent characteristics. This model is used to 
accurately simulate the effect of process variation and environmental changes on bias- 
dependent characteristics. The data- fitting model is used to model these characteristics 
with relatively good fidelity. The expressions of the data- fitting model are constructed 
to be charge conservative. As such, the resultant model is physically representative of 
the device is computationally robust within the harmonic balance algorithms employed 
by known large signal microwave circuit CAD tools. 

DESCRIPTION OF THE DRAWINGS 

These and other advantages of the present invention will be readily understood 
with reference to the following specification and attached drawings wherein: 

FIG. 1 is schematic diagram of an exemplary small signal equivalent circuit 
model for a HEMT device. 

FIG. 2 is a sectional view of an exemplary HEMT illustrating the rough 
translation of the physical origins for each of the equivalent circuit elements illustrated 
in the small signal circuit model in FIG. I. 

FIG. 3 is a cross-sectional view of a HEMT illustrating the regions in the HEMT 
which correspond to the various circuit elements in the small signal equivalent circuit 
model illustrated in FIG. I. 

FIG. 4 is an example of a cross-sectional descriptive input of a physical HEMT 
device structure by a conventional physical device simulation tool. 

FIG. 5 is an example of an epi stack descriptive input of a physical HEMT device 
structure for a known physical device simulation tool. 

FIG. 6 is an example illustrating the location of the epi stack within the device 
structures cross-sectional view. 

FIG. 7 is an example of a relatively accurate measured-to-model I-V 
characteristics using the semi-physical modeling method in accordance with the present 
invention. 

FIG. 8 is a elevational view illustrating an epi stack for an exemplary HEMT. 
FIG. 9 is a cross-sectional view of a HEMT for the exemplary epi stack illustrated 
in FIG. 8. 

FIG. 10 is a blown up diagram of the cross-sectional parameters pertaining to the 
T-gate geometry for the exemplary epi stack illustrated in FIG. 8. 
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FIG. 1 1 is a diagram of an electric conductance model used in the semi-physical 
example. 

FIG. 12 is a Smith chart illustrating the measured vs modeled S-parameters S 1 1, 
S12 and S22 simulated in accordance with the method in accordance with the present 
invention. 

FIG. 13 is similar to FIG. 12 and illustrates the measured vs. modeled values for 
the S21 parameter. 

FIG. 14 is similar to FIG. 12 but for the S12 S-parameter. 

FIG. 1 5 is a graphical illustration illustrating the fitting of non-linear data-fitting 
model to a semi-physical ly modeled I-V characteristics in accordance with the present 
invention. 

FIG. 16 is a graphical illustration illustrating an example of the measured-to- 
model I-V characteristics of a semi-physical modeling method. 

FIG. 1 7 is an implementation of the nonlinear data- fitting model as a user defined 
model in LIRBA 6.1. 

FIGS. 1 8 A and 1 8B illustrate the measured vs. modeled input power as a function 
of output power for a 22-23 GHz high-power amplifier. 

FIG. 1 9 is a block diagram i Uustrating the semiconductor semi-physical modeling 
in accordance with the present invention. 

FIG. 20 is a block diagram of the semi-physical model used in the processed 
perturbation to measured model method for modeling semiconductors which utilizes 
S-parameter microscopy in accordance with the present invention. 

FIG. 2 1 A is a schematic cross-sectional diagram of a standard HEMT used in the 
PM 2 experiment. 

FIG. 21 B is a cross-sectional diagram illustrating the epistack for the exemplary 
HEMT device used to demonstrate the model. 

FIG. 22 is schematic diagram illustrating the correspondence of the small signal 
equivalent circuit components to the detail of the device physical structure. 

FIG. 23 is a schematic diagram of the source access conductance of the HEMT. 

FIG. 24A is a graphical illustration of a source resistance I^ as a function of the 
biasing voltage V gs for different drain-to-source voltages V ds . 
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FIG. 24B is a graphical illustration of the source resistance R. as a function of the 
gate-to-source voltage illustrating the measured vs. semi-physically modeled 
approaches. 

FIG. 25 represents an exemplary S-parameter microscope in accordance with the 
present invention. 

FIG 26 illustrates the internal and external regions of an exemplary HEMT 
device. 

FIG. 27 is similar to FIG. 25 but illustrates the approximate locations of the 
model elements in the HEMT FET device illustrated is FIG. 25. 

FIG. 28 is a schematic diagram of a common source FET equivalent circuit 
model. 

FIG. 29 is an illustration of specific application of the S-parameter microscope 
illustrated in FIG. 25. 

FIG. 30 is similar to FIG. 25 which demonstrates the inability of known systems 
to accurately predict the internal charge and electrical field structure of a semiconductor 
device. 

FIG. 3 1 is a plan view of a four-fingered, 200|im GaAs HEMT device. 

FIG. 32 is a graphical illustration illustrating the measured drain-to-source current 
I Js as a function of drain-to-source voltage Vds for the sample FET device illustrated 
in FIG. 31. 

FIG. 33 is a graphical illustration illustrating the drain-to-source current I ds and 
transconductance G m as a function of the gate-to-source voltage V gs of the sample FET 
device illustrated in FIG. 31. 

FIG. 34 is a Smith chart illustrating the measured SI 1, S 12 and S22 parameters 
from frequencies of 0.05 to 40.0 GHZ for the FET device illustrated in FIG. 3 1 . 

FIG. 35 is a graphical illustration of the magnitude as a function of angle for the 
S21 S-parameter for frequencies of 0.05 to 40 GHz for the exemplary FET illustrated 
in FIG. 31. 

FIG. 36 is a graphical illustration of a charge control map of the charge and 
electric field distribution in the on mesa source access region shown with as a 
function bias in accordance with the present invention. 
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FIG. 37 is a graphical illustration of a charge control map of charge and electric 
field distribution in the on-mesa drain access region shown with R,, as a function of bias 
in accordance with the present invention. 

FIG. 38 is a graphical illustration of a charge control map for the non-quasi static 
majority carrier transport, shown with as a function of bias in accordance with the 
present invention. 

FIG. 39 is a graphical illustration of a charge control map for gate modulated 
charge and distribution under the gate, shown with Cgs and Cgt as function of bias in 
accordance with the present invention. 

FIG. 40 is a plan view of an exemplary Tt-FET with two gate fingers. 

FIG. 41 is a plan view of a Tt-FET with four gate fingers. 

FIG. 42 is an illustration of a tt-FET parasitic model in accordance with the 
present invention. 

FIG. 43 is an illustration of an off-mesa parasitic model for a Tt-FET in 
accordance with the present invention. 

FIG. 44 is an illustration of an interconnect and boundary parasitic model in 
accordance with the present invention for the Tt-FET with four gate fingers as illustrated 
in FIG. 41. 

FIG. 45 is an illustration of an inter-electrode parasitic model in accordance with 
the present invention. 

FIG. 46 is a schematic diagram of the inter-electrode parasitic model illustrated 
in FIG. 45. 

FIG. 47 is an illustration of an on-mesa parasitic model in accordance with the 
present invention. 

FIG. 48 is a schematic diagram of the on-mesa parasitic model illustrated in FIG. 

47. 

FIG. 49 is an illustration of an intrinsic model in accordance with the present 
invention. 

FIG. 50 is a schematic diagram of the intrinsic model illustrated in FIG. 49. 
FIG. 51 A is an exemplary device layout of a tt-FET with four gate fingers. 
FIG. 5 IB is an equivalent circuit model for the Tt-FET illustrated in FIG. 37A. 
FIG. 52 is a single finger unit device cell intrinsic model in accordance with the 
present invention. 
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FIG. 53 is similar to FIG. 52 and illustrates the first level of embedding in 
* accordance with the present invention. 

FIG. 54 is similar to FIG. 52 and illustrates the second level of embedding in 
accordance with the present invention. 
5 FIG. 55 is an equivalent circuit model of the 7t-FET illustrated in FIG. 51 A in 

accordance with the present invention. 

FIG. 56 is similar to FIG. 54 and illustrates the third level of embedding in 
accordance with the present invention. 

FIG. 57 is similar to FIG. 54 and illustrates the fourth level of embedding in 
10 accordance with the present invention. 

FIG. 58 is similar to FIG. 54 and illustrates the fifth level of embedding in 
accordance with the present invention. 

FIG. 59 A and 59B is a flow chart of a parameter extraction modeling algorithm 
that forms a part of the present invention. 
1 5 FIG. 60 and 6 1 illustrate the evaluation of a global error metric involving filtering 

error and trial model solution convergence speed. 

FIG 62A is a Smith chart illustrating the measured versus the initial model 
solutions for the SI 1, S 12 and S22 S-parameters from frequencies from 0.05 to 40.0 
GHz. 

20 FIG. 62B is a graphical illustration of angle versus magnitude for the initially 

modeled S-parameter S21 from frequencies of 0.05 to 40GHz. 

FIG. 63A is a Smith chart illustrating the measured versus simulated S- 
parameters SI I, SI 2 and S22 for frequencies 0.05 to 40 GHz for the first extraction 
optimization cycle. 

25 FIG. 63B is a graphical illustration of magnitude as a function of angle for the 

measure and first optimized model S-21 parameter for frequencies 0.05 to 40 GHz for 
the first optimization cycle. 

FIG. 64A is a Smith chart illustrating the measure as a function of the final model 
solution for S-parameters S 1 1 , S 1 2 and S22 for frequencies 0.05 to 40 GHz for the final 
30 solution. 

FIG. 64B is a graphical illustrations of the magnitude as a function of an angle 
for S-parameter S21 for the final model solution from frequency 0.05 to 40 Ghz. 
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FIG. 65 is a graphical illustration of the semi-physically modeled vs measured 
small signal Gm. 

FIG. 66 is a graphical illustration of the semi-physically simulated bias 
dependence of the small-signal output conductance Rds. 

FIG. 67 is a graphical illustration of the semi-physically simulated bias 
dependence of the small signal gate-source and gate-drain capacitance Cgs and Cgd. 

FIG. 68 is a graphical illustration of the semi-physically simulated bias- 
dependence of the small signal gate source charging resistance Ri. 

FIG. 69 is a graphical illustration of the semi-physcially bias dependence of the 
small signal source and drain resistance Rs and Rd. 

FIG. 70 is a graphical illustration of the measured vs modeled bias dependence 
gain at 23.5 Ghz for a K-band MMIC amplifier. 

FIG. 71 A and 71 B are graphical illustrations of the extracted parameters from 
measured device I-V's for process control monitor testing. 

FIG. 72 is a graphical illustration of the measured vs semi-physically simulated 
process variation for Gmpk and Vgspk. 

FIG. 73 is a graphical illustration of the measured vs semi-physically simulated 
process variation for Idpk and Gmpk. 

FIG. 74 is a graphical illustration of the measured vs semi-physically simulated 
process variation for Imax and Vpo. 

FIG. 75 is a graphical illustration of the measured/extracted vs semi-physically 
simulated process variation for the small signal equivalent model Rds and Gm. 

FIG. 76 is a graphical illustration of the measured/extracted vs semi-physically 
simulated process variation for the small signal equivalent model Cgs and Gm. 

FIG. 77 is a graphical illustration of the measured vs semi-physically simulated 
physical dependence for Imax as a function of physical gate length. 

FIG. 78 is a graphical illustration of the measured/extracted model vs semi- 
physically simulated physical dependence for Rds as a function of physical recess 
undercut width. 



DETAILED DESCRIPTION 

The present invention relates to a hybrid model based upon using a semi-physical 
model in conjunction with a data-fitting nonlinear model for use in large signal 
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microwave circuit computer aided design (CAD). The goal of the combined or hybrid 
modeling approach is to combine the physical modeling capability of the semi-physical 
approach along with the ability to be used in harmonic balance-based circuit simulation. 
The basic approach is set forth below: 

1 . Derive a semi-physical device model that is able to replicate measured DC I-V 
and bias-dependence small signal characteristics very accurately. 

2. Select a data-fitting nonlinear device model whose expressions satisfy the 
charge conservation law. 

3. Use a semi-physical device model to generate the necessary I-V 
(current/voltage) in C-V (capacitance/voltage) data for the data-fitting nonlinear device 
model. 

4. Derive model parameters for the data-fitting model such that it fits the semi- 
physical ly modeled characteristics. 

5. Implement the data-fitting device model in a large signal microwave circuit 
CAD tool. 

Steps 3-5 may be repeated for arbitrary physical changes to the semi-physical 
device. The semi-physical device model is able to reflect the effect of these changes 
in terms of the I-V and C-V performance. Consequently, through the use of data-fitting 
model, this change in performance can be simulated in a large signal regime. 

As mentioned above, Step 1 of the procedure is to derive a semi-physical device 
model for the HEMT device. The semi-physical device model is discussed below in 
connection through FIGS. 1-14. The hybrid model is discussed in connection with 
FIGS. 15-18. 

SEMI-PHYSICAL DETERMINATION OF SMALL-SIGNAL EQUIVALENT 

CIRCUITS 

The semi-physical model provides model elements for the standard small signal 
equivalent circuit model or FET as illustrated in FIG. 1 . However, unlike conventional 
methods, the model elements are derived from small signal excitation analysis of the 
intrinsic charge and electric fields within the device. As such, the simulated small 
signal model elements represent a relatively accurate physical equivalent circuit 
description of a physical FET. 
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The general methodology for the semi-physical modeling of intrinsic charge, 
electrical conductance and electric field is as set forth below: First, the relationships 
between the conduction band offsets, electric permitivity and material composition for 
the various materials in the epi stack are determined. These relationships can be 
performed analytically or by fitting simulated data from physical simulators. 
Subsequently, the basic electron transport characteristics in any of the applicable bulk 
materials in the epi stack are determined. Once the electron transport characteristics are 
determined, the undepleted linear channel mobility is determined either through 
material characterization or physical simulation. Subsequently, the Schottky barrier 
height value or expressions are determined. Once the Schottky barrier height value is 
determined, the semi-physical equations are constructed modeling the following 
characteristics: 

Fundamental-charge control physics for sheet charge in the active channel as 
controlled by the gate terminal voltage. 

Average centroid position of the sheet charge within the active channel width. 

Position of charge partitioning boundaries as a function of gate, drain and source 
terminal voltages. 

Bias dependence of linear channel mobility in the surface depleted regions. 
Bias dependent of the velocity saturating electric field. 
Saturated electron velocity. 

Electrical conductance of the linear conductance region of the channel under the 

gate. 

Electrical conductance within the source and drain access regions 
Once the semi-physical equations are determined, the empirical terms of the semi- 
physical modeling equations are adjusted to fit the model I-V (current/voltage) 
characteristics against measured values. Subsequently, the empirical terms are 
iteratively adjusted to achieve a simultaneous fit of measured C-V (capacitance-voltage) 
and I-V characteristics. Lastly, the empirical modeling terms are fixed for future use. 

By constructing a comprehensive set of semi-physical equations that cover all of 
the physical phenomenon as mentioned above, the physical operating mechanisms 
within a HEMT device can be relatively accurately determined. FIG. 7 illustrates a set 
of relatively accurate measured-to-modeled I-V characteristics for a HEMT using the 
semi-physical modeling discussed herein. In particular, FIG. 7 illustrates the drain-to- 
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source current 1^ as a function of the drain-to-source voltage for various gate biases, 
for example, from 0.4V to -1 .0V. As shown in FIG. 7, solid lines are used to represent 
the semi-physical model while the Xs are used to represent measured values. As shown 
in FIG. 7, a close relationship exists between the measured values and the modeled 
5 parameters. 

An example of semi-physical modeling for physical device operation in 
accordance with the present invention is provided below. The example utilizes an 
exemplary device as illustrated in FIGS. 8 and 9. Table 2 represents exemplary values 
for the physical cross-section dimension parameters of the device. FIG. 1 0 relates to 
10 a blown up T-gate characteristic which is correlated to the parameters identified in 

Table 1. 
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Table 1 

Values for the Physical Parameters Input into Device Cross Section 



Layout Parameter 


Units 


Value 


Gate Length 




[Mm] 


0,150 


Wing Length 


Lgw 


[Mm] 


0.520 


Gate Mushroom Crown Length 


Lgmcl 


[Mm] 


0.200 


Total Gate Height 


Hg 


[Mm] 


0.650 


Gate Stem Height 


Hgstem 


(wnl 


0.300 


Gate Sag Height 


Hgsag 


[Mm] 


0.100 


Gate Cross-Sectional Area 


Gate Area 


[Mm 2 ] 


0,187 


Max Cross-Sectional Area 


MaxArea 


[Mm 2 ] 


0.364 


Total Gate Periph 


Wg 


[Mm] 


200.000 


# Fingers 


N 


[1 


4.000 


Source-Drain Spacing 


Osd 


[Mm] 


1.800 




Dsg 


[Mm] 


0.700 


Gate-Drain Spacing 


uga 


fMm] 


1.100 


Gate-Source Recess 


RECsg 


[Mmj 


0.160 


Gate-Drain Recess 


RECgd 


[Mm] 


0.240 


Recess Etch Depth 


Hrec 


[A] 


780.000 


SiN Thickness 


Hsin 


[A] 


750.000 


Gatefeed-Mesa Spacing 


Dgfm 


[Mm] 


2.000 


Gateend-Mesa Overlap 


Dgem 


[Mm] 


2.000 


Finger-Finger Spacing Thru Drain 


Dffd 


[Mm] 


16.500 


Finger-Finger Spacing Thru Source 


Dffs 


[Mm) 


13.500 


Source Airbridge Inset? 


AB? 


[] 


P 


Source Airbridge Inset 


Dsabin 


[Mm] 


28.000 


Source Airbridge Height 


Hsab 


[Mm] 


3.500 


Source-Gate Airbridge Clearance 


Hgsab 


[Mm] 


1.640 


Source Pad Width 


Ws 


[Mm] 


12.000 


Drain Pad Width 


Wd 


[Mm] 


14.000 


Substrate Thickness 


Hsub 


[Mm] 


100.000 



As mentioned above, the semi-physical modeling of the intrinsic charge and 
5 electric field within the HEMT device is initiated by determining the relationships 

between the conduction band offset, electric permitivity, and material composition for 
the various materials in the epi stack. Material composition related band offset and 
electrical permitivity relationships may be obtained from various references, such as 
"Physics of Semiconductor Devices," by Michael Shur, Prentice Hall, Englevvood 

10 Cliffs, New Jersey 1990. The basic electron transport characteristics, for example, for 

the l inear mobility of electron-carriers in the bulk GaAs cap layer may be determined 
to be 1350cm 2 /Vs, available from "Physics of Semiconductor Devices", supra. The 
linear mobility of electron carriers in the undepleted channels is assumed to be 
5 500cm W s. This value may be measured by Hall effect samples which have epi stacks 

1 5 grown identically to the stack in the example, except for some differences in the GaAs 

cap layer. The Schottky barrier height is assumed to be 1 .051 volts, which is typical of 
platinum metal on a AlGaAs material. 
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HEMT device. These expressions govern how to partition the model charge between 
the influence of different terminals. 



Effective Gato Lengtn 


L BOtf 




= L 0 *aL.* A L« 


Gpitp.Soitrrfl Pnntmf Rornon 

vjoic*uuun,t3 vUl III W V%ISQPUil 




\vM 




Source-Side Effective G3le Length Extension 




U m l 


= A»-.e * ' V 


Orain-Side Effective Gate Length Extension 




l(i m l 




Gate-Drain Control Region 




1^1 


= (L^* i L - )*{«anh|10{LV2.X 0 ,J|*l}/2 








■ XotV^M^/ 








{M xrfK V at /(U[X 0L V,.M x<JL / 


Bias Dependent Extension of me Saturated Transport Region 








Empirical Ora in-Saturated Transport Boundary Factor 






= <- g V i( /{2|lMV^./n 








= 1<(1»(VJV. - jnr*». 








Vd.«"<v - ,/v, Mn r ,| /{ 








V.„ n m(t + (V <1 ./V.., n n ,,,m *"}1 


Position of the Boundary between Regions 1 and 2 


x s 






Norto: Region t denotes ma 1 in eer region, wn 


lie R«tjlo« 2 denote* (tie saturated region of the channel 


Empirical Effective Gate Length Extension Gate Bias Factor 




t M m/V] 




Empmcat Effective Gate Length Extension Drain Bias Factor 




(um/Vl 




Effective Dram-Source VoltageCcntrot-2 








Rough, fntrinsic Saturation Vottage 




M 












Rough, intrinsic Saturation CuoeM Level 


'sal 


[A| 


[i * 0c*R, ♦ sqn( Wg^'(V Bt ,A/;n) 


intrinsic Conductance of the Lnear Region. Under the gate 


9ch> 


(SJ 


S W N. ^.Wg)/Lg 


Rough Intrinsic Saturation voiage Level 


Vl 


m 




Empmcat Knee Shaping Parameter 


m 


[] 




Empirical Region 2 extension Orain Bias Factor 


M Xdt 


(i 




Empincal Reg:on 2 extension Gate Bias Factor 


M XdK 


(] 




Fine intnnsic Saturation Voltage 


V s»t»n 




= 1 /g ^ 








= g cw v flt . v L "(-v t (A*a*R.)* 








sqrt{V l i {A*g 6f ,,R,)^V Btt J .(g wl R 1 V l } I )| / 


Fine uunnsic Saturation CujtwiI Level 


' vi team 


[A] 


fv„,. a < '-g c ^{v L /v 01 .)-,i 


Saturation Region Length Ratio 


A 


() 


= Xs/L,., 


initial Starting position for Region 1 & 2 Boundary 




l M m ) 


= Lg/2 


Region 1 & 2 3 oundary Bos Factor 


Mx. 


(1 




Region i & 2 Boundary Drain Bos Factor 
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Region 1 £ 2 Boundary Drain Bias Factor 


MxsK 


(] 





The following equations represent the semi-physical expressions used to model 
the bias dependence of linear channel mobility in depleted regions. 

Depleted Channel Mobility tL3V9 {cm 2 ATs] = Mdehan ♦ m -k -V„. 

Fixed Depleted Channel Mobflfty ^i dC han l crTl2/v *sJ 
Depleted Channel Mobility Gate Bias Factor ^ dK lcm 2 /V**s] 



The following equations are the semi-physical expressions used to model the bias 
dependence of the saturating electric field and saturation velocity. 



15 



WO 01/84466 



PCT/USO 1/13361 



Saturating Electric Field 
Fixed Saturating Channel Mobility 
Saturating Channel Mobility Gate Bias Factor 
Saturation Velocity 



Usat 
UsatK 



[V/um] 
[cm 2 /V*s] 
[cm*/V 2 *s] 

[cm/sj 



= v, / 
K^sat + ^Vgt.) 10000J 



FIG. 1 1 is a schematically illustrates how electrical conductance in the source and 
drain access regions are modeled in the example. 

The following equations describe the semi-physical model for the source access 
region conductance. 



Source Access Resistance 




tnJ 


Source Access Rp<;i%(ancp- Channel »«rl <~ar\ 


o 

S undapCap 




Source Access Resistance: Recess and Undeleted Cap 


R S4cc»u 




Source Access Resistance. Crowding resistance due to conductance mismatch 


^ SBoundwry 


tn. M m] 


Resistance of the Source Recess Access region at high cn-state bias (V^) 


a ON 
n S d*pR«c 


(0. M m) 


Resistance of the Undeleted Source Recess Access region 


p 

S undopRtc 




Uncapped. Fully Depleted Sheet Resistance 


^SHd«p 


(n/sq) 


Capped. Undepleted Sheet Resistance 


R SH 


[n/sql 


Uncapped. Unoepieted Sheet Resistance 


^SHurtdvp 


Cn/sq] 


Cap Sheet Resistance 


R SHC*p 


In/sqi 


Surface Oepietion Factor 


^sunVndvp 


[] 


High On-state bias: Diode Turn-on voltage 


V 0N 


tv) 


Ohmic Contact Resistance 


R cont 




RF Ohmic Contact Resistance Reduction Factor- 


RF rconF 


[] 


Source Access Resistance Bias Modification Factor 


MR, 


U 


Cf-Vds Bias Modification Factor 


MC,, 


[J 


Rs-Vds Bias Modification Factor 


MR.L 


[] 


Rs-Vgs Bias Modification Factor 


MR. K 


U 


Cf-Vgs Swith point to On-state 


vc^ 


(I 


Cf-Vgs Bias Expansion Factor 


KC m 


[] 


Rs-Vgs Swith point to On-state 


VRson 


[) 


Rs-Vgs Swith point to Off-state 


VRsorr 


(1 


Rs-Vds Swith pcmt from Off-On transition 


VR SKn«. 


[] 


Rs-Vds Bias Expansion Factor 


XR.L 


(] 


Rs-Vgs. Bias Expansion Factor 


KR, K 


11 


Rs-Vds Bias Expansion Factor @ Rs Saturation 


KR iSal 


[I 


Rs Bias Shaping Factor 




[] 



~ ( Rsinrf«vC^ * RsAocm« * RsBoundary ) ' ^fl 
= Re^RF (eonf ♦ 

'WO aa -<REC lf *L l /2H 
= R&^~'MR/tanh<[KC m - ( 
V«. - VC^V^'MC^ JJ ♦ 1 ) / 2 • { V fl . 

/ 2 -p-tanh< KR, K (V 0 . - VR.oj)!)* ( tanh { 

KR, Srt (V - ,-VR„ u# .)j*i )/2 ♦ R SvntfM ^ • { lanht 

XR.k(V 0 .-VR. O(1 )I*1)/2 
= R S d^ ON 'MR,*ianh(|KR. K -V^*KR. t -V d .+ 
VR.«r I ♦ 1 ) / 2 ' ((1 ♦V^MR, J'MR, K *j 1 - 

tanh(KR. ni O/ < ,-VR, K(1# .))J/ | 2'(1 ♦ [V„./«1 ♦ 

s Rsh^CREC^) 

= RsKmdap ( R EC^ ) 

= 1 M 1/R S mc*p * t/ R SMund . p ) 



= F.. 



= t/ lQ N^, M CM) ( H c-p -H c , 



The following equations describe the drain access region conductance: 
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Drain Access Resistance 




[Qj 


Drain Access Resistance, except for recess access 


"ORF 


InJ 


vjratn rvecfSas access r^estsiancs 


Q 


lol 


Drain Access Resistance: Recess and Undepteted Cap 


" D*ee«ss 




urain Access rvestsiance. ^n«nnet ana v_*ap 


R O undepCap 


U L w }L ■ I 


Kesistancs or utb unoepieieo urain r\etess Access region 




rn«nm1 


Resistance of the Saturated Drain Recess Access region 


" DSsturatit 




Resistance of the Drain Recess Access region at high on-state bias (V„) 


a on 
"OdepR»c 






"ofl«* 




Drain Access Resistance Bias Modification Factor 


MR d 


(] 


Rd-Vds Access Bias Modification Factor 


MR daccL 


[] 


Rd'Vgs Swith point to On -state 


v **ttOn 


[] 


Rd-Vgs Swith point at saturation 


VR dSat 


CI 


Rdaccess-Vgs Bias Expansion Factor 


KRdaecK 


U 


Rd-Vds Saturation Bias Modification Factor 


MR dSat 


tl 


Rd-Vds Access Resistance Bias Modification Factor 


MR dL 


[] 


Rd-Vgs Saturation Swith point to Off-state 


WdSaiOff 


n 


Rd-Vgs Saturation Swith point 


VRdSat 


n 


Rd-Vgs Saturation Swith point transition width 


VR dSato 


tj 


Rd Bias Shaping Factor 


yRfl 


[i 



* ( Roup * Ron ) 
* Ron 

= Ro*«,fl^ ON (MR/v 0 .- 
tanh< »«W [ V,. - VR W „ ♦ V d 
: R«.«yRF rewrf ♦ 

= RsKuotf^ ( REC 9 „ 



.•MRd^.Ll)W2 

Rsh( D fld 

) 



= MR aSat %.exp{ -I V 8 ,W„-MR dL *VR dS . t0ff j 2 / (2 
VR d5rtw ) J / { VR -S „ [1 ♦< V dB /VR dJ- )' R - J< * } 
= RsMd«» ( REC gd ) 
s 5vl M ♦ < *W... ' Rtw.pc. P V( 26vt) * SOrtf 6 2 
♦ ((Rou... - Ro^pc«)/(26vt) - 1ft ) 



SEMI-PHYSICAL DETERMINATION OF 
SMALL-SIGNAL EQUIVALENT CIRCUITS 

To derive values for the familiar small signal equivalent circuit as shown in FIG. 
1, a small signal excitation analysis must be applied to the semi-physically modeled 
physical expressions. The method of applying such an analysis is as follows: 

1) Gate Terminal Voltage Excitation 

Apply a small*/- voltage delta around the desired bias condition, across the gate- 
source terminals. 

Equivalent circuit element Gm - delta(Ids)/delta (Vgs') where delta (Vgs*) is 
mostly the applied voltage deltas, but also subtracting out that voltage which is dropped 
across the gate source access region, shown as RsCont, RsundepCap, RsundepRec, 
ResdepRec, and RsBoundary in FIG. 11, above. 

Equivalent circuit element Cgs and Cgd takes the form of 
delta(Nsn)/de!ta(Vgs)*Lgn, where delta (Nsn) is the appropriate charge control 
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expression, and Lgn is the gate source or gate drain charge partitioning boundary 
length. 

Equivalent circuit element Ri=Lgs/(Cgschannel * vs) where Cgs channel is the 
portion of gate source capacitance attributed to the channel only, and vs is the saturated 
electron velocity. 

2) Drain Terminal Voltage Excitation 

Apply a small +/- voltage delta around te same bias condition as in 1 , but the delta 
is applied across drain source terminals. 

Equivalent circuit element Rds = l/{delta(Ids)/delta(Vds , )} where Vds 1 is mostly 
the applied voltage deltas, but also subtracting out voltage which is dropped over both 
the gate source and gate drain access regions. 

Equivalent circuit element Cds is taken to be the sum of the appropriate fringing 
capacitance Semi-Physical models, or can take the formofdelta(Nsd)/delta(Vds')*Xsat, 
were Nsd is the charge control expression for charge accumulation between the 
appropriate source and drain charge boundaries, and Xsat is the length of the saturated 
region, if in saturation. 

3) On-mesa Parasitic Elements: The equivalent circuit elements, Rs and Rd are 
expressed by the appropriate electrical conduction models of the source and drain 
access regions. 

The RF performance can be predicted at an arbitrary bias point. 

Table 3 represents a comparison of the values for a high frequency equivalent 
circuit model derived from equivalent circuit model extraction from and semi-physical 
modeling for the sample illustrated in Table 2. 



Table 2 

Comparison of Modeled Equivalent Circuit Results for Semi-physical 
Modeling Method, and Equivalent Circuit Model Extraction 



Intrinsic Equivalent 
Circuit Parameter 


Equivalent Circuit 
Model 


Semi-Physical 
Devic Model 


Cgs 


0.227745 pF 


0.182 pF 


Rgs 


64242 Q, 


infinite Q 
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Intrinsic Equivalent 


Equivalent Circuit j 


Semi-Physical 


Circuit Parameter 


Model | 


Device Model 


Cgd 


0.017019 pF 


0.020 pF 


Rgd 


133450 Q 


infinite Q 


Cds 


0.047544 pF 


0.033 pF 


Rds 


160.1791 Q 


178.1 a 


Gm 


135.7568 mS 


124 mS j 


Ri 


3.034 Q 


2.553 Q f 


Tau 


0.443867 pS 


0.33 pS 



The results of the semi -physical modeling method produce a small-signal 
equivalent circuit values which are relatively more accurate than the physical device 
simulator in this case. Furthermore, given the differences in the parasitic embedding, 
treatment of the two approaches, the results given in Table 3 yield much closer results 
than a comparison of equivalent circuit values. 

Table 3 lists the values of parasitic elements used in the model derivations. An 
important difference between the extracted equivalent circuit model and the semi- 
physically derived one is the use of Cpg and Cpd to model the effect of launch 
capacitance for the tested structure. This difference leads to the results of the extracted 
model results being slightly off from the optimum physically significant solution. 
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Table 3 

Comparison of Modeled "Parasitic" Equivalent Circuit Results f r Semi 
physical Modeling Method, and Equivalent Circuit Model Extraction 



Extrinsic Equivalent 


I Equivalent I 


I Semi-Physical 


circuit K3ramet©r 


circuit Mooei 


uevice Moaei 


Rg 


1.678 Q 


1.7Q 


Lg 


0.029314 nH 


0.03 nH 


Rs 


1.7Q 


1.21 Q 


Ls 


0.002104 nH 


0.003 nH 


Rd | 


3.309899 Q 


5.07 Q 


Ld | 


0.031671 nH 


0.02 nH 


Cpg | 


OpF 


0.02 pF 


Cpd 


OpF 


0.01 pF 
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As shown in FIGS. 12, 13 and 14, the modeled results that are simulated using 
the semi-physically derived equivalent circuit model very accurately replicate the 
measured high frequency, S-parameter data. 

The following equations represent the small-signal excitation derivation of 
20 small-signal equivalent circuit modeled Gm. FIG. 65 illustrates the semi-physically 

simulated bias equations of the small signal Gm compared to measured data. 



Semi- Physically Modeled Drain-Source Current Control l tf , {A| P + (V 4U /V Mu r)'"" 

Small-Signal Determination of equiv Circuit Gm value 9 m pF (S] s / d( V^. - v^, ) 



Source-Access voltage drop Vjacc M e 1 

Fine Extrinsic Saturation Voltage V irt% (Vj = t^^/g.* 

extrinsic Conductance of the Linear Region. Under tne gate g eh {SI (1 ♦g,*, (R,* R «)1 

Intrinsic Conductance of the linear Region .Under trie gate Q eM [5] =(Q N. 



The following equations represent the small-signal excitation derivation of Rds. 
FIG. 66 illustrates the semi-physically simulated biase-dependence of the small-signal 
25 Rds. 
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Small-Signal Deterimination of equiv. Circuit Rds value R ds fol = 1 / 9 d *RF 

= { d l d , / d[ V d , - R pfob#0 J d . - V BAce -V dA «-V d!Ut ) 
> *W 

= *d«W * ( R 0Und«pC«p + R 0Und*pR«c + R 0Acc«s» + 
Rprcb-O^Q ) 
= 'd»W * ( RoS*tor»t»d ) 







9d6RF 


[S] 


VpAcc 


[V) 


V 0SAT 


IV] 


RprobaD 


(nl 


»dsF 


[] 


R^rdsF 


[I 



= (RFrd»F + D* 

ianh( 10*|V d .-V lh |) + 1 



The following equations may be used for the small-signal excitation derivation 
of Cgs and Cgd. FIG. 67 illustrates the semi-physically simulated bias-dependence of 
the small-signal Cgs and Cgd. 



Small-Signal Deterimination of equiv. Circuit Rds value R ds 

9dsRF 

Drain-Access voltage drop V OAcc 

Drain-Saturated Region voltage drop Vqsat 

External Test probe or lead resistance R pr0 beD 

High Frequency conductance dispersion factor r dsF 

High Frequency conductance dispersion RP r dsF 



= { * l d , / d( V d , - R probi0 . I dl - V fAec -V dAee -V dSat ) 
IS} > * r dtF 

= *d*W * ( ^DUndopCjp + RoUnd«pR«c + ^Otecifi + 
(V] Rprob.0^0 ) 

[V] = * dtV V * ( RoS«tof«t#d ) 

[Qj 

= (RFrd.F + D* 

[] tanh( 10-|V d ,-V,J)+ 1 

[1 



The following equations are involved in the small-signal excitation derivation 
of RI FIG. 68, which follows, shows the semi-physically simulated bias-dependence 
of the small-signal Ri. 



Gate-Source Non-quasistatic charging resistance Rjchargo fo*M m l = L gs2 W g / [ C osChan W fl v s ] 

Gate-Channel Source Capacitance C giChan [fF/ M mJ { t . r (V Mton -V dM ) / ( 2*V, Jt . n -V d5e )] 2 } 



FIG. 69 below shows the semi-physically simulated bias-dependence of the on- 
mesa parasitic access resistances, Rs and Rd. 
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EXAMPLE OF SEMI-PHYSICAL MODEL AND BIAS-DEPENDENCE 
SMALL-SIGNAL SOURCE AND DRAIN RESISTANCE. RS AND RP. 

The following example verifies how the semi-physical small-signal device 
5 model is able to provide accurate projections for bias-dependent small-signal 

performance. In this example, the same semi-physical device model as used in the 
previous examples was used because the example MMIC circuit was fabricated utilizing 
the same HEMT device technology. 

In this example, the bias-dependence small-signal gain and noise performance 
of a two-stage balanced K-band MMIC LNA amplifier is replicated through microwave 
circuit simulation using small signal and noise equivalent circuits that were generated 
by the semi-physical model. The results of the measured and modeled results are 
shown below in Table 4. As seen from these results, the semi-physical device model 
was able to accurately simulate the measured bias-dependent performance, even though 
the bias variation was quite wide. 



Table 4 

Measured vs. Modeled Gain NF and Gain @ 23.5 Ghz for K-band MMIC LNA at 

Difference Bias Conditions 



Bias Condition 


Measured Gain 


Predicted Gain 


Measured NF 


Predicted 




@ 23.5 GHz 


@ 23.5 GHz 


@ 23.5 GHz 


NF @ 23.5 
GHz 


Vds=0.5 V 


15.2 dB 


15.8 dB 


2.97 dB 


2.77 dB 


112 mA/mm 










Vds=1.0 V 


20.6 dB 


21.0 dB 


2.29 dB 


2.20 dB 


112 mA/mm 










Vds=2.0 V 


19.8 dB 


20.2 dB 


2.25 dB 


2.15 dB 


112 mA/mm 










Vds=3;0 V 


18t9 dB 


19.1 dB 


2.30 dB 


2.11 dB 


112 mA/mm 










Vds=3.5 V 


18.4 dB 


18.5 dB 


2.34 dB 


2.18 dB 


112 mA/mm 










Vds=4.0 V 


18.0 dB 


18.0 dB 


2.37 dB 


2.27 dB 


112 mA/mm 










Vds=2.0 V 


16.4 dB 


18.0 dB 


2.45 dB 


2.21 dB 
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56 mA/mm 










Vds-2.0 V 


21.4 dB 


20.9 dB 


2.38 dB 


2.21 dB 


170 mA/mm 










Vds=2.0 V 


22.2 dB 


21.0 dB 


2.65 dB 


2.6 dB 


225 mA/mm 










Vds=3.0 V 


21.4 dB 


20.3 dB 


2.71 dB 


2.61 dB 


225 mA/mm 










Vds=3.0 V 


20.5 dB 


20.0 dB 


2.42 dB 


2.22 dB 


170 mA/mm 










Vds=4.0 V 


19.6 dB 


19.2 dB 


2.50 dB 


2.29 dB 


170 mA/mm 











A plot of measured vs. modeled gain for the values listed in Table 3, above, is 
shown in FIG. 70. 

EXAMPLE OF SEMI-PHYSICAL MODEL ACCURATE PHYSICAL 

NATURE 

The following example verifies how the semi-physical small-signal device 
model is able to provide accurate projections for physically dependent small-signal 
performance. Tn this example, the same semi-physical device model as used in the 
previous examples was used. 

In this example, physical process variation was input into the semi-physical 
device model in terms of statistical variation about known averages, cross-correlation, 
and standard deviations. The goal of this exercise was to replicate the measured DC 
and small-signal device variation. The degree of accurate replication indicates the 
degree to which the semi-physical model is physically accurate. 

Table 5 below lists the simulated, and known process variation that was used: 
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Table 5 

Statistical Process Variation Model 



Parameter 


Nominal 


Standard Dev. 


Gate Length 


0.15 urn 


0.01 um 


Gate-Source Recess 


0.16 um 


0.015 um 


Gate-Drain Recess 


0.24 um 


0.020 um 


Etch Depth 


780 A 


25 A 


Pass. Nitride Thickness 


750 A 


25 A 


Gate-Source Spacing 


0.7 um 


0.1 um 


Source-Drain Spacing 


1.8 um 


0.15 um 



In the course of microelectronic component production, sample devices are tested 
in process in order to gain statistical process control monitor (PCM) data. FIGs 71 A and 
7 1 B shows schematically the kind of data that is extracted and recorded from measured 
device I-V's during PCM testing. 

Since the semi-physical device model is able to simulate I-V's, it was able to 
simulate the variation of I-V's due to physical process variation. These I-V's were 
analyzed in the same fashion to extract the same parameters that are recorded for PCM 
testing. Figures 72, 73 and 74 show how accurately the simulated results match with 
measured process variation. Figure 72 shows how the semi-physically simulated Vgpk 
and Gmpk match with actual production measurements. Figure 73 shows how 
simulated Idpk and Gmpk match, also. Finally, FIG. 74 shows how simulated Imax and 
Vpo also match very well. 

Small-signal S-parameter measurements are also taken in process for process 
control monitoring. These measurements are used to extract simple equivalent circuit 
models that fit the measured S-parameters. Since the semi-physical device model is 
able to simulate these equivalent circuit models, it was able to simulate the variation of 
model parameters due to physical process variation. 

Figures 75 and 76 show how accurately the simulated results match with 
measured/extracted process variation for the small-signal model parameters. Figure 75 
shows how the semi-physically simulated Rds and Gm match very well with actual 
extracted model process variation 
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More direct and convincing evidence supporting the accurate, physical nature of 
the semi-physical model can be shown be comparing the dependence of simulated and 
measured performance to real physical variable. As shown in Figure 77, the semi- 
physical model is able to very accurately reproduce the dependence of Imax upon gate 
length. In addition, the semi-physical model is also able to replicate physical 
dependence for high-frequency small-signal equivalent circuits. This is shown in 
Figure 78, which shows that it is able to reproduce the dependence of Rds with Recess 
undercut width. 

HYBRID APPROACH 

As mentioned above, the modeling approach in accordance with the present 
invention is a hybrid of semi-physical and data-fitting modeling. The first step in the 
hybrid modeling approach is to derive a semi-physical device model as discussed 
above. Secondly, a data-fitting nonlinear device model whose expressions satisfy the 
charge conservation law is selected. An exemplary data-fitting nonlinear device model 
is disclosed in "A New Empirical Model for HEMT Devices/' by Y. C. Chen, IEEE 
Microwave Guided Wave Letters . Vol. 8, No. 10, October 1998. This model satisfies 
the charge conservation constraint for current. In order to complete the description of 
the nonlinear model the following charge expressions are used: 

=(C 8inax +C fimin ) * V ss / 2 + C ed * (V^-VJ 
Gate Charge (Empirical Nonlinear Model) Q G [C] (C BttW - C fimin ) * In {cosh [a c (V gs -V g0 )] } / (2\) 



This equation satisfies the conservation constraints for charge as follows: 

dVgsdVds dVdsdVgs 
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This condition, which is requisite for path integration independence of the 
expressed charge is satisfied by the expression above. 

FIG. 15 demonstrates that the chosen nonlinear data-fitting model relatively 
accurately fits the semi-physically modeled I-V characteristics. The nonlinear model 
that is chosen in this example is "A New Empirical Model for HEMT Devices, supra. 
FIG. 16 demonstrates how the semi-physical model relatively accurately models the 
measured IV characteristics. By transitivity the nonlinear data-fitting model fits the 
measured characteristic. 

Table 6 below lists the model parameters used with the data- fitting model to fit 
the semi-physically modeled data. 
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Table 6 

Parameters Used in the Nonlinear Data-Fitting Model 



Data-Fitting 
Model 

pa ram pfpr 

rai ai i icici 


Units 


Value 


Wref 


[um] 


! 200 


Tau 

--— -y^ ------- ----- 

Cgmax 

"CgmTn 

------- — — — — — — — — 

Acg 

-------••■J----—-—— 

Vcg 
cds 
CcTg 
Crf 


[PS] 

— -fpFj" 


| 0.35 
0.178 

0:03 

"575 


TV] 

— "1 fpFJ 
CpJT" 


I -ff.73 

Q[irT6Z67 

TflOOOW 


Rc 


TohmJ 


5000 


Kl 


[onm] 


5.841 fS3873 


Ffg 


TohmJ 


0.858T33333 


Ffs 


ToKml 


1338f2S967 


RaT 


[onm J 


317323603? 


Lg 


TnH] 


| 0.BFTB 


L"d 


TnHI ] 


CUTS' 


Ls 


TnH] 


0 


Is aid i ode 


TA] 


TE-T3 ] 


Ncflocf© 




f.2 


RcfTocfe 


"foKmJ 


32V5 


Tpfc 


TA] 


TE 1 22875" 


VR 


L ^ 


TE6CS207 


Lam 




(F.0f9G5'S 


ADO 


u 


^3. 167863 


*ADf 




£088357 


AD2 


U | 


4T00f5U7 


A10 






mi 


B 1 


^J. (3 1*9706 




B | 


^J.0053~6? 


A20 


U J 


-T.2032TJ? i 




H J 


57167313 


1 


H 




-AW 


Tl 

U 


— rarfi^ — 


"A3T 


— F — 




"A32 


-Q- 


~"TnJ50SW"~ 


7«U 


"""U — 


""""03r48"655""" 


TJUPf 


— "XT— ' 


— qCT3B5S5*~ 


AW 


F — 


— &OT8B3T"" 




--TV]"" 


"076 


v gmin 


""TV] — 


-rz — 


" Tflagf * 


"Tohmr 


100OOW"" 


■ — -mr g *s" — 


"TohmJ" 


35000 


Vpo 


"TV1""~ 


— "3:6797""" 



FIG. 17 illustrates the implementation of the nonlinear data- fitting model into a 
commercially available microwave CAD tool, for example, LIBRA 6.1 by Agilent 
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Technologies. The LIBRA 6. 1 CAD tool contains large signal simulation capabilities 
based on the harmonic balance algorithm. The particular view shown in FIG. 17 
illustrates the circuit element schematic that was implemented and not the actual 
enabling implementation of the user defined model which is normally implemented as 
a computer program code extension to LIBRA 6. 1 . The procedures by which user- 
defined models may be added to LIBRA 6.1 are discussed in detail in User's Manual 
entitled "User-Defined Elements", Series VI User Manuals, released 6.0 Hewlett 
Packard, HP Part No. E4605-90033, July 1995. 

FIGS. 1 8A and 1 8B illustrate how the hybrid device modeling produces relatively 
accurate large signal simulation results. In particular, the input power response is 
shown for a two-stage balanced high-power amplifier at a frequency of 22.5 GHz. The 
measured responses are shown in FIG. 18B while the simulated response is shown in 
FIG. 18 A. As shown, the measured vs. model responses are relatively close to each 
other indicating a relatively high degree of accuracy for the device modeling process. 

The following example demonstrates the ability of the approach to accurately 
incorporate the effect known physical process changes into the simulation of large 
signal device characteristics. The multiple curves that are shown in each plot of FIGS. 
1 8 A and 1 8B may be obtained from: several measured circuits from a plurality (i.e. 1 3) 
wafers in the case of measured data; and simulated performance variation due to 
simulated physical process variation. The simulated process variation used is as 
discussed below under the caption "PM 2 Process Perturbation to Measured Model 
Method for Semiconductor Device Technology Modeling". The parameters used to 
simulate the process variation are provided below in Table 7. 

The tables that have been generated for a 0.15 micrometer HEMT process relate 
to a realistic distribution of physical process variation. The following process 
parameters were varied accordingly as shown in Table 7. 
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Table 7 

Physical Process Variation Model Parameters 



Parameter 
Gate Length 
Gate-Source Recess 
Gate-Drain Recess 
Etch Depth 

Pass. Nitride Thickness 
Gate-Source Spacing 
Source-Drain Spacing 



0.15 urn 
0.16 um 
0.24 um 
780 A 
750 A 
0.7 um 
1.8 um 



Nominal 



25 A 
0.1 um 
0.15 um 



0.01 um 
0.015 um 
0.020 um 
25 A 



Standard Dev. 



A physical correlation may be simulated between the recess spacing and gate 
length. A correlation factor of 0.5 between gate length and both recess dimensions may 
be incorporated in the semi-physical model statistics. Although the process variation 
described above is simple and realistic it is adequate in describing all of the DC and RF 
performance variations for wafer means/medians. 

PROCESS PERTURBATION TO MEASURED MODELED METHOD FOR 
SEMICONDUCTOR DEVICE TECHNOLOGY MODELING 

A method for modeling semiconductor devices based upon a process perturbation 
to measured modeled (PM 2 ) methodology can be used to develop a physically-based 
technology model that ultimately becomes more and more accurate as more and more 
process perturbation experiments are performed. As shown in FIG. 19, various 
parameters, such as device scaling, bias dependence, temperature dependence, layout 
dependence and process dependence can be modeled using this technique to analyze 
measurements taken for any imaginable set of process perturbations. The more 
measurements that are taken, the semi-physical model becomes more and more 
"corrected". For example, by performing more PM 2 experiments in which the gate 
length of high electron mobility transistor (HEMT) samples are varied to much longer 
lengths than originally studied, the models for velocity saturation and effective gate 
source charge control length can be refined to provide more accurate results for longer 
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gate lengths. Also, by performing temperature dependent measurements the 
temperature dependence on the material parameters is able to be refined to better fit the 
modeled to the measured results. 

An important part of the PM 2 modeling methodology is a measured-to-model 
microscope (i.e. S-parameter microscope) which is able to look into the "guts" of a 
semiconductor device. With this capability a relatively comprehensive physically- 
based model for the entire device technology can be developed. 

The perturbation method is discussed below in connection with FIGS. 19-24. An 
important aspect of the perturbation method is a measured-to-model microscope (S- 
parameter microscope) as discussed below in connection with FIGS. 25-39. The 
measured-to-model microscope may utilize a filter in order to remove the contribution 
of device layout parasitics to the modeled electrical characteristics. This may be done 
to accomplish clear representations of the internal physical operation for the measured 
device. One embodiment of such a filter for P-FET layouts is discussed in connection 
with FIGS. 40-58. The exemplary measured-to-model microscope, also called an S- 
parameter microscope, utilizes an extraction algorithm for extracting modeled 
parameters as generally discussed in connection with FIGS. 59-64. 

The following example illustrates the use of the PM 2 modeling concept to create 
a complete, physically-based model for the source resistance of a HEMT device. The 
PM 2 experiment is used to determine the physical model characteristics are as follows: 

1 ) Characterize standard HEMT samples 

A) Use standard fabrication process to produce HEMT samples that are embedded 
within a standard device layout. 

B) Collect information regarding the physical dimensions of the source-access 
region by way of a scanning electron microscope (SEM). 

C) Test the sample devices using S-parameter microscopy to establish physically 
representative equivalent models. 

2) Characterize the standard HEMT samples with device layout experiments 
A) Use standard fabrication process to produce HEMT samples that have device 

layout experiments that vary the physical dimensions of the source access region, for 
example, gate source spacing etc. 
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B) Collect information regarding the physical dimensions of the source access 
region by way of SEM. 

C) Test the sample devices with S-parameter microscopy to establish physically 
representative equivalent models. 

3) Characterize the HEMT samples with a thin GaAs "cap" 

A) Use standard fabrication process to produce thin cap HEMT samples 
embedded within the standard device layout. 

B) Collect information regarding the physical dimensions of the source access 
region by way of a SEM. 

C) Test using S-parameter microscopy to establish physically representative 
equivalent models. 

A cross-sectional diagram of the HEMT sample used in the example is illustrated 
in FIG. 21 A. A diagram of a standard device for a standard Pi-FET example layout is 
shown in FIG. 5 1 A. A cross-sectional diagram of the material epi-stack present in the 
HEMT is shown in FIG. 7B. For the third part of the PM 2 experiment above, the GaAs 
cap is thinned down to 7.5 nm instead of the standard thickness of 50 nm keeping the 
same doping density. 

After a sample wafer is fabricated for the standard HEMT sample, Scanning 
Electron Microscopy (SEM) is used to determine the dimensions of the critical 
structural components. The measured and intended structural dimensions are identified 
in Table 8 below, where measured refers to SEM determined dimension and standard 
indicates the normal or intended specification. Each of the dimensions listed in Table 
8 are correlated to the cross-sectional diagram illustrated in FIG. 21 A. 
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Table 8 

SEM Measured and Intended Structural Dimension for a "Standard"HEMT 

Device Sample 





units 


Measured 


Standard 


Lg 


[Mm] 


0.15 


0.15 


Dsg 


[umj 


0.78 


0.8 


Dsd 


[N 


1.935 


2 


Lg+RECsg+RECgd 


[Mm] 


0.56 


0.52 



Next, S-parameter microscopy as described below and in connection with FIGS. 
25-64 is used to determine physically representative, model representations of the 
source access resistance. The exemplary Pi-FET may be modeled and used as a filter 
in the S-parameter microscopy as generally described. This is accomplished by 
measuring the S-parameters of the sample devices up to 40 GHz and subsequently 
extracting equivalent small signal circuit models as discussed in detail below in 
connection with S-parameter microscopy. The small signal equivalent circuit model 
serves as an electrical representation of the physical structure of the measured device 
and can be used to roughly sketch the details of its internal structure. The 
correspondence between equivalent circuit elements and structural items within the 
device are shown in FIG. 28 below. The relationship of the quantity "Rs" and the 
source access region is shown. 

The results of the S-parameter microscopy measurements are shown in FIG. 36 
which plot the bias dependent characteristics of the source resistance Rs. From these 
bias dependent characteristics, a preliminary physical model which fits the measured 
data can be constructed. 

Three physical effects were found to contribute to the overall behavior of the 
measured source resistance: resistance of the access region before the region; resistance 
within the source access recess; and a boundary resistance caused by sudden change in 
sheet carrier concentration between the source access and the channel directly under the 
gate. These phenomena and their physical locations within the source access recess are 
illustrated in FIG. 23, where the regions 1, 2 and 3 correspond to each of the effects 
discussed above. 
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From these observations, a form of semi-physical model describing bias- 
dependent characteristics of Rs in a HEMT device can be established as set forth in the 
equations below. 



1 S undepCap 



- ( RsundepCap + RsAccess + RsBoundary ) I W g 

= Rcon/R^rconF + 
fo*^ m J Rsh( - (REC sg +Lg/2)] 



In part two of the PM 2 experiment, the HEMT device samples are fabricated and 
tested and the length of the source access region is intentionally varied. After the 
samples are fabricated, the intended dimensions are verified through SEM. S-parameter 
microscopy is also used to extract the source resistance for comparison. The 
experimental source gate dimensions along with the extracted source resistance at a bias 
point of Vds = 1 .0, Vgs = 0.01 are provided in Table 9 below. 

Table 9 

Extracted Source Resistance vs. Source-Gate Distance 



[Dsg ^m] ||Rs [q] 



Measured Sheet Res: 
Extracted Sheet Res: 



0.4 

0.S 


1 

o 


1 

r:i 


TV*5 


109.1 
110 


n/Sq -1 
hTsq 



The data in Table 9 is used to confirm the preliminary semi-physical model for 
Region 1 's source-access resistance (Rsundep Cap) illustrated above. This confirmation 
can be verified by comparing the extracted sheet resistance (Rsh) by S-parameter 
microscopy and the PM 2 experiments against sheet resistance extracted by an 
independent Van der Pauw measurements. "Modem GaAs Processing Methods"" by 
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Ralph Williams, Artech House, Norwood, MA 1 990. Even though the experiment may 
be conducted using HEMT devices with a different material or epi stack, the experiment 
illustrates the validation of a semi-physical model form for Region 1 resistance. Also 
the terms RECsg and Lg may be assumed to be roughly constant for all of the Dsg test 
samples. In the final part of the PM 2 experiment, the full form of the semi-physical 
source resistance model is validated. Based on the full bias dependent measurement of 
the part 1 , the complete semi-physical model expression source resistance as a function 
of gate and drain bias can be represented in the equations below: 



Rs [qI ~ ( ^SundepCap + RsAccess + RsBoundav ) ' 

KsuncfepCap ^Q*^ - (REC^2)] 

= Rs^Rec^^MRs^tanhflKC^* ( 

V gs / 2 11-tanh( KR, K (V gs - VR^))]}* { 

Access k>^] { tanh[ KRj K (V^-VR^)] + 1} / 2 

= ^ep^^-M^^tanh^K^K^H.KR,^^ 
VR^ ] + 1 } / 2 * {(HV^J^k^I 

te"H^Ms^Knee)W I 2*0 + [^((1 

SBoundary 



The simulated result for the sample fabricated in part I of the PM 2 experiments 
is shown in FIG. 24A. Comparing FIG. 24A with FIG. 36 indicates that the semi- 
physical model adequately replicates the measured results. The result of the experiment 
is shown below in FIG. 24B. As expected the bias dependent source resistance of the 
thin "cap" sample has the same form, only offset higher by an amount that corresponds 
to the change in Rsh in Region 1 of the source access. 
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S-PARAMETER MICROSCOPY 

The S-parameter Microscopy (SPM) method utilizes bias dependent S-parameter 
measurements as a form of microscopy to provide qualitative analysis of the internal 
charge and electrical field structure of the semiconductor device heretofore unknown. 
Pseudo images are gathered in the form of S-parameter measurements extracted as 
small signal models to form charge control maps. Although finite element device 
simulations have heretofore been used to calculate the internal charge/electric field of 
semiconductor devices, such methods are known to be relatively inaccurate. The S- 
parameter microscopy provides a relatively accurate method for determining the 
internal charge and electric field within a semiconductor device. With accurate 
modeling of the internal charge and electric field, all of the external electrical 
characteristics of the semiconductor devices can be relatively accurately modeled 
including its high frequency performance. Thus, the system is suitable for making 
device technology models that enable high frequency MMIC yield analysis forecasting 
and design for manufacturing analysis. 

S-parameter microscopy is similar to other microscopy techniques in that SPM 
utilizes measurements of energy reflected to and from a sample to derive information. 
More particularly, SPM is based on transmitted and reflective microwave and 
millimeter wave electromagnetic power or S-parameters. As such, S-parameter 
microscopy is analogous to the combined operation of scanning and transmission 
electron microscopes (SEM and TEM). Scattered RF energy is analogous to the 
reflection and transmission of the electron beams in the SEM and TEMs. However, 
instead of using electron detectors as in the SEM and TEMs, reflectometers in a 
network analyzer are used in S-parameter microscopy to measure a signal. S-parameter 
microscopy is similar to other microscopy techniques in that both utilize; measurement 
of scattering phenomenon as data; include mechanisms to focus measurements for 
better resolution; and include mechanisms to contrast portions of the measurement to 
discriminate detail as shown in Table 10 below: 
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Table 10 



General Microscopes 


S-Parameter Microscope 


Measure of scattered energy 


Measures S-Parameters 


Mechanism for "focus" 


Focuses bv extraction of Unique 
equivalent circuit models 


Mechanism for "contrast" 


Contrasts by using bias dependence to 
finely discriminate the nature and 
location of charge/electric fields 



RESULT : Detailed "images" of device's internal charge and electric field 
structure. 

Images as discussed herein, in connection with S-parameter microscopy do not 
relate to real images, but are used provide insight and qualitative detail regarding the 
internal operation of a device. More specifically, S-parameter microscopy does not 
provide visual images as in the case of traditional foims of microscopy. Rather, S- 
parameter microscopy images are more like maps which are computed and based on a 
non-intuitive set of measurements. 

FIG. 25 illustrates a conceptual representation of an S-parameter microscope 
generally identified with the reference numeral 20. The S-parameter microscope 20 is 
analogous to a microscope which combines the principles of SEM and TEM. Whereas 
SEM measures reflections and TEM measures transmissions, the 2-port S-parameter 
microscope 20 measures both reflective and transmitted power. As a result, data 
derived from the 2-port S-parameter microscope contains information about the 
intrinsic and extrinsic charge structure of a device. More particularly, as in known in 
the art, SEM provides relatively detailed images of the surface of a sample through 
reflected electrons while TEM provides images of the internal structure through 
transmitted electrons. The reflective signals are used to form the external details of a 
sample while transmitted electrons provide information about the interior structure of 
a device. In accordance with an important aspect of the invention, S-parameter 
microscopy utilizes a process of measuring reflective and transmitted signals to provide 
similar "images" of the charge structure of a semiconductor device. As used herein the 
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internal and external electrical structure of a semiconductor device are commonly 
referred to as intrinsic device region and 22 and extrinsic parasitic access region 24 as 
shown in FIG. 26. Also contributing to the external electrical structure of the device 
are parasitic components associated with its electrodes and interconnects, which are not 
shown. These are "layout parasitics" the so-called device or its so-called internal 
electrical structure. 

Referring to FIG. 25, the ports 26 and 28 are emulated by S-parameter 
measurements. The S-parameter measurements for a specific semiconductor device, 
generally identified with the reference number 30, are processed in accordance with the 
present invention to provide charge control maps, shown within the ci rcle 32, analogous 
to images in other microscopy techniques. These charge control maps 32, as will be 
discussed in more detail below, are expressed in the form of equivalent circuit models. 
As shown in FIG. 27, the linear circuit elements may be used in the models to represent 
the magnitude and state of charge/electric fields inside the semiconductor device 30. 
The position of the circuit elements within the model topology is roughly approximate 
the physical location within the device structure, hence the charge control map 
represents a diagram of the device's internal electrical structure. 

The interpretation of the exact location of measured charge/electric fields within 
the semiconductor device is known to be ambiguous since an equivalent circuit model, 
for example, as illustrated in FIG. 28 with discrete linear elements, is used to represent 
the distributed structure of the charge/electric fields in the actual device. Although 
there is no exact method for distinguishing the physical boundaries between measured 
quantities, bias dependence is used to clarify how the S-parameters should be 
discriminated, separated and contrasted. In particular, changing bias conditions is 
known to change the magnitude and shift boundaries between the charge and electric 
fields within the device. The changes are normally predictable and qualitatively well 
known in most technologies. As such, the charge control maps can readily be used as 
maps illustrating the characterization of physical changes in magnitude, location and 
separation of electric charge and electric fields. 

Analogous to other forms of microscopy, the S-parameter microscope 20 
emulates a lens, identified with the reference numeral 40 (FIG. 25). The lens 40 is 
simulated by a method for the extraction of a unique equivalent circuit model that also 
accurately simulates the measured S-parameter. More particularly, parameter extraction 
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methods for equivalent circuit models that simulate S-parameters are relatively well 
known. However, when the only goal is accurately fitting measuring S-parameters, an 
infinite number of solutions exist for possible equivalent circuit parameter values. 
Thus, in accordance with an important aspect of the present invention, only a single 
unique solution is extracted which accurately describes the physical charge control map 
of the device. This method for unique extraction of equivalent circuit model parameters 
acts as a lens to focus the charge control map solution. As discussed and illustrated 
herein, the lens 40 is subsequently simulated by a filter that is based on an apparent 
layout parasitic embedding model. As discussed below, the layout parasitic embedding 
model consists of linear elements which simulate the effect of the device's electrodes 
and interconnects upon its external electrical characteristics. A Pi FET embedding 
model 42, is described below. This model effectively acts as a filter to remove the 
electrical structure of the extrinsic parasitic access contribution to the preliminary 
charge control map solution. The resultant filtered charge control map solution 
represents a clearer "image" which shows only the electrical structure of the intrinsic 
device. This enhanced imaging is needed in order to achieve as accurate a view of the 
internal electric charge/field as possible. Unlike conventional extraction techniques as 
illustrated in FIG. 30, which can only extract non-unique equivalent circuit models and 
not the unique charge control map, the S-parameter microscope 20 is able to relatively 
accurately model the internal electric charge/field structure within a semiconductor 
device. 

An exemplary application of the S-parameter microscope is illustrated in detail 
below. In this example, an exemplary GaAs HEMT device with four gate fingers and 
200|am total gate periphery formed in a Pi-FET layout as generally illustrated in FIG. 
31 and identified with the reference numeral 43, is used. The GaAs HEMT 43 is 
adapted to be embedded in a 1 00-fim pitch coplanar test structure to facilitate on water 
S-parameter measurement. 

Initially, as illustrated in FIGS. 32 and 33, the I-V characteristics for the device 
are measured. In particular, the drain source current Ids is plotted as a function of 
drain-to-source voltage Vds at various gate voltages Vgs as shown in FIG. 32. and FIG. 
33 illustrates the drain-to-source current Ids as a function of gate voltage Vgs and 
transconductance Gm (i.e. the derivative of Ids with respect to Vgs) at different drain 
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voltages Vds. These I-V characteristics are typical of HEMT devices, which are one 
type of three-terminal semi-conductor device technology. 

Table 1 1 shows the bias conditions in which S-parameters were measured. The 
S-parameters were measured from 0.05 to 40 GHz at each bias condition. FIG. 34 
illustrates a Smith chart illustrating the measured S-parameters SI 1, S 12 and S22 for 
frequencies from 0.05 to 40.0 GHz. FIG. 35 is a graphical illustration of magnitude as 
a function of angles for the measured S-parameter S2 1 for frequencies from 0.05 to 40.0 
GHz. 

TABLE 11 
Measured S-parameter Bias Conditions 



Biases 

Vgs 


Vds = 
0 V 


Vds = 
0.5 V 


Vds = 
1.0 V 


Vds = 
2.0 V 


Vds = 
4.0 V 


Vds = 
5.0 V 


-1.6 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


-1.4 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


-1.2 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


-1 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


-0.8 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


-0.6 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


-0.4 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


-0.2 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


0 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


0.2 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


0.4 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


0.6 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 1 



Using the small signal model illustrated in FIG. 28, the extracted small signal 
equivalent circuit values are obtained as illustrated in Table 12 for each S-parameter at 
each bias condition, using the extraction method discussed below. 
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The values in Table 1 2 represent solutions that are close to the charge control map 
and represent physically significant solutions of the FET's electrical structure. 
However, the values represented in Table 12 contain the influence of external layout 
parasitics which, in accordance with an important aspect of the invention, are subtracted 
using a model for the embedding parasitics to obtain the most accurate charge control 
mapping to the intrinsic device characteristic. In particular, an embedding model is 
applied to filter the extracted equivalent circuit model values and obtain values more 
representative of the intrinsic device. In particular, in the exemplary embodiment, a 
PiFET embedding parasitic model is used to subtract capacitive contributions due to 
interelectrode and off-mesa layout parasitic influences. This filter essentially subtracts 
known quantities formed from the parameters Cgs, Cgd and Cds depending on the 
device layout involved. In this example, embedding of the inductive parameters is not 
necessary because these quantities are extrinsic and do not contribute to the charge 
control map of the intrinsic device. 

As discussed above, the lens with filter is used to generate unique charge control 
maps. In particular, FIGS. 36-39 illustrates the bias dependent charge control maps for 
the parameters RS, RD, RI, CGS and CGD as a function of bias. More particularly, 
FIG. 36 illustrates a charge control map of the charge and electric field distribution in 
the on-mesa source access region illustrated by the source resistance as a function 
of bias. FIG. 37 illustrates a charge control map of the charge and electric field 
distribution in the on-mesa drain access region illustrated by the drain resistance R d as 
a function of bias. FIG. 38 illustrates a charge control map for a non-quasistatic 
majority carrier transport illustrated by the intrinsic device charging resistance Rj as a 
function of gate bias for different drain bias points. FIG. 39 illustrates a charge control 
map for gate modulated charge and distribution under the gate shown with the gate 
capacitance CGS and CGD as a function of bias. 
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FILTER 

As mentioned above, the S-parameter microscope 20 utilizes a filter to provide 
a clearer charge control map for modeling the internal electric charge/field of a 
semiconductor device. Although the filter is illustrated in connection with the PiFET 
with multiple gate fingers, as illustrated in FIGS. 40 and 41, the principles of the 
invention are applicable to other semiconductor devices. 

As illustrated in FIG. 40, PiFETs are devices in which the gate fingers and the 
edge of the active region resemble the greek letter ic, as illustrated. Such PiFET layouts 
facilitate construction of multi fingered large periphery device cells, for example, as 
illustrated in FIG 41. In accordance with an important aspect of the invention, the 
multi-finger semiconductor device is modeled as a combination of single finger device 
cells. Each single finger device cell is represented by a hierarchy of four models, 
which, in turn, are assembled together using models for interconnects to represent an 
arbritrary multi-fingered device cell, illustrated in FIG. 42. The four models are as 
follows: off mesa or boundary parasitic model; interelectrode parasitic model; on-mesa 
parasitic model and intrinsic model. 

The off-mesa parasitic model is illustrated in FIG. 43. This model represents the 
parasitics that exist outside the active FET region for each gate finger. In this model, 
the fringing capacitance of each gate finger off the active device region as well as the 
off-mesa gate finger resistance is modeled. 

The interelectrode parasitic model and corresponding equivalent circuit are 
illustrated in FIGS. 44-46. This model represents parasitics between the metal 
electrodes along each gate finger. The following fringing capacitance parasitics are 
modeled for the gate-to-source air bridge; drain-to-source air bridge; gate-to-source 
ohmic; gate-to-drain ohmic and source-to-drain ohmic as generally illustrated in FIG. 
45. 

The on-mesa parasitic model and corresponding equivalent circuit are illustrated 
in FIGS. 47 and 48. This model represents that parasitics around the active FET region 
along each gate finger including various capacitance fringing parasitics and resistive 
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parasitics. In particular, the gate-to-source side recess; gate-drain-side recess; gate- 
source access charge/doped cap; and gate-drain access charge/doped cap capacitance 
fringing parasitics are modeled. In addition, the gate metallization and ohmic contact 
resistive parasitics are modeled. 

The intrinsic model and corresponding equivalent circuit are illustrated in FIGS. 
49 and 50. The intrinsic model represents the physics that predominately determine the 
FET performance. In particular, the DC and current voltage response can be 
determined by physics based analytical equations for magnitude and location of 
intrinsic charge which are generally know in the art, for example, as disclosed in 
"Nonlinear Charge Control in AlGaAs/GaAs Modulation-Doped FETs," by Hughes et 
al, IEEE Trans. Electron Devices, Vol. ED-34, No. 8, August 1987. The small signal 
model performance is modeled by taking a derivative of the appropriate charge or 
current control equations to derive various terms such as RJ, RJ, RDS, RGS, RGD, GM, 
TAU, CGS, CDS and CGD. Such control equations are generally known in the art and 
disclosed in detail in the Hughes et al reference mentioned above, hereby incorporated 
by reference. The noise performance may be modeled by current or voltage 
perturbation analysis "Noise Characteristics of Gallium Arsenride Field-Effect 
Transistors" by H. Statz, et al IEEE-Trans. Electron Devices , vol. ED-21, No. 9, 
September 1974 and "Gate Noise in Field Effect Transistors at Moderately High 
Frequencies" by A. Van Der Ziel, Proc. IEEE , vol 51, March 1963 and "Noise 
Characteristicsof Gallium Arsenide Field Effect Transistors," by H. Statz, IEEE Trans. , 
Vol. ED-21, No. 9, September 1974. 

Ail example of a parasitic model for use with the S-parameter microscopy 
discussed above is illustrated in FIGS. 51-58. Although a specific embodiment of a 
semiconductor device is illustrated and described, the principles of the present invention 
are applicable to various semiconductor devices. Referring to FIG. 51 A, a Pi-FET is 
illustrated. As shown, the PiFET has four gate fingers. The four fingered Pi-FET is 
modeled in FIG. 5 IB. In particular, FIG 5 IB illustrates an equivalent circuit model for 
Pi-FET illustrated in FIG. 5 1 A as implemented by a known CAD program, for example, 
LIBRA 6. 1 as manufactured by Agilent Technologies. As shown, the equivalent circuit 
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models does not illustrate all of the equivalent circuit elements or network connections 
involved with implementing the parasitic embedding models, but rather demonstrates 
a finished product. FIG. 5 IB is displayed in symbol view in order demonstrate 
resemblance to FIG. 41. The actual technical information regarding the construction 
of the network and its equivalent circuit elements are normally provided in schematic 
view. 

FIGS. 52-56 demonstrate the application of the parasitic model for use with the 
S-parameter microscopy. An important aspect of parasitic modeling of multi-gate 
fingered devices as single gate finger devices. As used herein, a single unit device cell 
refers to a device associated with a single gate finger. For example, a four fingered Pi- 
FET as illustrated in FIG. 51 is modeled as four unit device cells. 

Initially, the four finger Pi-FET illustrated in FIG. 51, is modeled as a single 
finger unit device cell 100 with an intrinsic model 102, as shown in FIGS. 52 and 53. 
In particular, the Pi-FET intrinsic FET model 104 is substituted for the block 102 
defining a first level of embedding. As shown in FIG. 53, the parameter values for the 
Pi-FET intrinsic model are added together with the parameter values for the single 
fingered unit device cell intrinsic model. The intrinsic device model 104 may be 
developed by S-parameter microscopy as discussed above. Next, as illustrated in FIG. 
54, the interconnect layout parasitic elements are added to the equivalent model by 
simply adding the model terms to the value of the appropriate circuit element to form 
a single unit device cell defining a second level of embedding. Once the single unit 
device cell is formulated, this device is used to construct models for multi- fingered 
devices. In this case, a Pi-FET with four gate fingers is modeled as four single finger 
device unit cells as shown in FIG. 55. Subsequently, the off-mesa layout parasitic 
elements are connected to the multi-fingered layout, defining a third level of embedding 
as illustrated in FIG. 56. These off-mesa layout parasitic elements, generally identified 
with the reference numerals 108 and 110, are implemented as new circuit elements 
connected at key outer nodes of the equivalent circuit structure. Subsequently, a fourth 
level of embedding is implemented as generally illustrated in FIG. 57. In particular, an 
inductor model is connected to the sources of each of the various unit device cells to 
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represent the metallic bridge interconnection, as generally shown in FIG. 54. Lastly, 
as illustrated in FIG. 58, a fifth level of embedding is implemented in which the feed 
electrodes model 1 14 and 1 16 are modeled as lumped linear elements (i.e. capacitors 
inductors) as well as the distributive elements (i.e. microstrip lines and junctions) to 
form the gate feed and drain connections illustrated in FIG. 58. As shown, the 
distributive elements are distributed models for microstrip elements as implemented in 
LIBRA 6.1. 

EXTRACTION METHOD FOR UNIQUE DETERMINATION OF FET 
EQUIVALENT CIRCUIT MODELS 

The method for determining FET equivalent circuit parameters as discussed 
above is illustrated in FIGS. 58-64. This method is based on an equivalent circuit 
model, such as the common source FET equivalent circuit model illustrated in FIG. 28. 
Referring to FIG. 58, a model is initially generated in step 122. The model illustrated 
in FIG. 28 is used as a small signal model for the FET. In accordance with an important 
aspect of the algorithm, the equivalent circuit parameters are based upon measured FET 
S-parameters. Measurement of S-parameters of semiconductor devices is well known 
in the art. FIG. 62A is a Smith chart illustrating exemplary measured S-parameters 
SI 1, SI 2 and S22 for frequencies between 0.05 to 40 GHz. FIG. 62B represents a 
magnitude angle chart for the measured S-parameter S21 from frequencies from 0.05 
to 40 GHz. After the S-parameters are measured, as set forth in step 124 (FIG. 59), it 
is ascertained whether the measurements are suitable in step 126. This is either done 
by manually inspecting the test result for anomalies, or by algorithms to validate the test 
set. If the measurements are suitable, the S-parameter measurements are stored in step 
128. 

A space of trial starting feedback impedance point values, for example, as 
illustrated in Table 13 is chosen. Then, a direct model attraction algorithm, known as 
the Minasian algorithm, is used to generate preliminary values for the equivalent circuit 
model parameters, for each value of starting feedback impedance. Such extraction 
algorithms are well known in the art, for example, as disclosed "Broadband 
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Determination of the FET Small Equivalent Small Signal Circuit" by M. Berroth, et al., 
IEEE - MTT, Vol. 38, No. 7, July 1990. Model parameter values are determined for 
each of the starting impedance point values illustrated in Table 14. In particular, 
referring to FIG. 59 A, each impedance point in Table 1 2 is processed by the blocks 1 30, 
132, etc. to develop model parameter values for each of the impedance point in order 
to develop an error metric, which, in turn, is used to develop a unique small signal 
device model, as will be discussed below. The processing in each of the blocks 130, 
132 is similar. Thus, only a single block 130 will be discussed for an exemplary 
impedance point illustrated in Table 13. In this example, the impedance point 1 7 which 
correlates to a source resistance R. ohm of 1 .7Q and a source inductance L s of 0.OO45pH 
is used. 

TABLE 13 

Trial Starting Feedback, Impedance Space Point Values 



Impedance 


1 Resistance 


Inductance 1 


Point 


1 (Rs) 


(Ls) 


1 


0.1 Q 


0.0045 pH 


1 2 


0.2 Q 


0.0045 pH 


i 3 


0.3 Q 


0.0045 pH 


— t — 


0.4 Q 


0.0045 pH 




0.5 Q 


0.0045 pH 


! 6 


0.6 Q 


0.0045 pH 


! 7 


0.7 Q 


0.0045 pH 


8 


0.8 Q 


0.0045 pH 


9 


0.9 Q 


0.0045 pH 


10 


1.0Q 


0.0045 pH 




1.1 Q 


0.0045 pH ! 


12 


1.2 Q 


0.0045 pH 


13 


1.3 n 


0.0045 pH 


14 


1.4 Q 


0.0045 pH 


15 


1.5 Q 


0.0045 pH 


16 


1.6 a 


0.0045 pH 


17 


1.7 Q 


0.0045 pH 


18 


1.8 Q 


0.0045 pH 


19 


1.9 Q 


0.0045 pH 


20 


2.0 Q 


0.0045 pH 
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21 


| 2.1 Q 


u.uuto [j n 


22 


2.2 Q 


n nr)4^ nM 


23 


2.3 Q. 


u.uu'w (Jn 


24 


2.4 a. 


0.0045 pH 


25 


2.5 a 


0.0045 pH 


26 


2.6 Q 


0.0045 pH 


f 27 


2.7 Q 


0.0045 pH 


28 


2.8 Q 


0.0045 pH 


29 


2.9 Q 


0.0045 pH 


30 


3.0 Q 


0.0045 pH 



For the selected value, R s = 1.7 ohms, initial intrinsic equivalent circuit 
parameters and initial parasitic equivalent circuit parameter are determined, for 
example, by the Minasian algorithm discussed above and illustrated in Tables 14 and 
15 as set forth in steps 134 and 136. In step 138 the simulated circuit parameters are 
compared with the measured S-parameters, for example, as illustrated in FIGS. 62. 
Each of the processing blocks 130 and 132 etc. goes through a fixed number of 
complete cycles, in this case six complete cycles. As such, the system determines in 
step 140 whether the six cycles are complete. 



TABLE 14 

Initial "Intrinsic" Equivalent Circuit Parameters 



Intrinsic Equivalent 


Initial Solution 


Circuit Parameter 




Cgs 


0.23595 pF 


Rgs 


91826 Q 


Cgd 


0.0177 pF 


Rgd 


100000 Q 


Cds 


0.04045 pF 


Rds 


142.66 Q 


Gm 


142. 1025 mS 


Tau 


0.1 pS 
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TABLE 15 

Initial "Parasitic" Equivalent Circuit Parameters 



Intrinsic Equivalent 




Circuit Parameter 


| Initial Solution 


Rg 


3.0 CI 


Lg 


0.014 nH 


Rs 


1.7 Q 


Ls 


0.0045 nH 


Rd 


2.5 Q 


Ld 


0.024 nH 



Each cycle of the processing block 1 30 consists of a direct extraction followed 
by an optimization with a fixed number of optimization iterations, for example 60. By 
fixing the number of extraction-optimization cycles along with the number of 
optimization iterations, a fixed "distance" or calculation time which the model solution 
must be derived is defined. As such, the algorithm implements a convergence speed 
requirement of the global error metric by setting up an environment where each trial 
model solution competes against each other by achieving the lowest fitting error over 
a fixed calculation time thus causing a "race" criteria to be implemented, where 
"convergence speed" to be implicitly calculated for each processing block 1 30, 1 32 etc. 

After the system determines whether the racing is done in step 140, the system 
proceeds to block 142 and optimizes model parameters. Various commercial software 
programs are available. For example, a commercially available, LIBRA 3.5 software 
as manufactured by HP-eesof may be used both for circuit simulation as well as 
optimizing functions. The optimization is performed in accordance with the restrictions 
set forth in Table 16 with the addition of fixing the feedback resistance R. to a fixed 
value. 
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TABLE 16 

Environment Used for Competitive Solution Strategy, as Implemented in this 

Example 



Implementation Parameter | 


Circuit Simulator and Optimizer 


Libra 3.5 


Optimization Algorithm 


Gradient 


Optimization Error Metric 


Mag and angle of 
Sll,S21,S12,and S22 
from 4 to 40 GHz 


Number of Iterations 


60 


Number of Extraction/Optimization 

I Cycles 


6 



By fixing the value for R. this segment of the algorithm confined to creating a 
trial model solution for only the trial feedback impendence point with which it started. 
Table 17 illustrates the optimized intrinsic equivalent parameter values using 
commercially available software, such as LIBRA 3.5. These values along with the 
optimized parasitic values, illustrated in Table 18, form the first optimized model 
solution for the first extraction-optimization cycle (i.e. one of six). The optimized 
model parameters are then fed back to the function block 134 and 136 (FIG. 63 A and 
63B) used for a new initial model solution. These values are compared with the 
measured S-parameter value as illustrated in FIG. 63. The system repeats this cycle for 
six cycles in a similar fashion as discussed above. After the six extraction-optimization 
cycle, the final trial model solution for the trial impendence point 17 is complete along 
with its final fitting error to the measured data to form the new error metric 144. In 
accordance with an important aspect, the extraction-optimization algorithm makes the 
final optimization fitting error for each point implicitly carry information about both the 
measured to model fitting error and the speed of convergence. It does so by the fixed 
optimization time constraint which sets up a competitive race between the various trial 
model solutions. 
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TABLE 17 

Optimized "Intrinsic 11 Equivalent Circuit Parameters 



I Intrinsic Equivalent 
Circuit Parameter 


1 Initial Solution 


Cgs 


0.227785 pF 


Rgs 


65247 Q 


Cgd 


0.017016 pF 


Rgd 


F 130820 Q 


Cds 


0.047521 pF 


Rds 


160.18 Q 


Gm 


135.74 mS 


Tau 


0.446 pS 


TABLE 18 

Optimized "Parasitic" Equivalent Circuit Parameters 


Intrinsic Equivalent 1 
Circuit Parameter 


Initial Solution 


Rg 


4.715 Q 


Lg 


0.02903 nH 


Rs* 


1.7 a 


Ls 


0.002102 nH 


Rd 


3.2893 Q 


Ld 


0.0317 nH 



The implementation of the extraction optimization cycles makes the best and 
fastest solving solution appear as a global minima for the final fitting error in step 146 
of all of the trial impedance points as generally shown in FIGS. 60 and 61. More 
specifically, referring to FIG. 60 the global minima solution using the new error metric 
is found around 1^=1.7 ohms. Tables 19 and 20 list the final model equivalent circuit 
parameters for this global solution, including the intrinsic and parasitic parameter as set 
forth in step 148 (FIG. 59). 
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TABLE 19 

Global Solution for "Instrinsic" Equivalent Circuit Parameters 



j Intrinsic Equivalent Circuit 


Initial Solution 1 


Parameter 




Cgs 


0.227745 pF 


Rgs 


64242 Ci 


Cgd 


0.017019 pF 


Rgd 


133450 Q 


Cds 


0.047544 pF 


Rds 


160.1791 n 


Gm 


135.7568 mS 


1 Tail 


0.443867 pS 



TABLE 20 

Global Solution "Parasitic" Equivalent Circuit Parameters 



1 Extrinsic Equivalent Circuit j 


1 Initial Solution 


j Parameter 




Rg 


4.711895 Q 


tg j 


0.029314 nH 


Rs 


1.7 O 


Ls | 


0.002104 nH 


Rd 1 


3.309899 Q 


Ld | 


0.031671 nH 



In order to test the accuracy of the solution, the final model for solutions are 
compared with the measured S-parameter values as shown in FIG. 64. As shown, there 
is good correlation between the simulated model values and the measured S-parameters 
values thus veri fy ing that the simulated model values represent a relatively accurate and 
unique small signal device model. 

Obviously, many modifications and variations of the present invention are 
possible in light of the above teachings. Thus, it is to be understood that, within the 
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scope of the appended claims, the invention may be practiced otherwise than as 
specifically described above. 

What is claimed and desired to be covered by a Letters Patent is as follows: 
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I CLAIM: 

1. A method for modeling semiconductor comprising the steps of: 

(a) developing a semi-physical model; 

(b) developing a data-fitting model; and 

(c) combining said semi-physical and data-fitting model to form a hybrid model. 

2. The method as recited in claim 1, wherein step (c) includes the following 

steps: 

(d) using the semi-physical model to generate current- voltage (I-V) and 
capacitance-voltage (C- V) data for the data fitting model; 

(e) deriving model parameters for the data-fitting model such that it fits the 
semi-physically modeled characteristics. 

3. The method as recited in claim 2, wherein step (c) further includes step (f): 
implementing the data-filling device model in a large signal microwave circuit 
computer aided drafting (CAD) tool. 

4. The method as recited in claim 3, further including step (g): repeating steps 
(d)-(f) for arbitrary physical changes to the semi-physical device. 

5. The method as recited in claim 1, wherein step (a) includes the step (h): (h) 
determining the relationships between the conduction band offsets, electric permitivity 
and material composition of the materials contained in the epi-stack. 

6. The method as recited in claim (5), wherein step (h) is performed 
analytically. 

7. The method as recited in claim (5) wherein step (h) is performed by fitting 
simulated data from physical simulators. 
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8. The method as recited in claim 5, wherein step (a) also includes step (i): 
determining the basic electron transport characteristics in the materials in the epi-stack. 

9. The method as recited in claim 8, wherein step (a) also includes step 
determining the undepleted linear channel mobility. 

10. The method as recited in claim 9, wherein step (j) is determined through 
material characterization. 

1 1 . The method as recited in claim 9, wherein step (j) is determined by physical 
simulation. 

12. The method as recited in claim 9, wherein step (a) also includes step (k): 
determining the Schottky barrier expressions. 

13. The method as recited in claim 1, wherein step (a) includes semi-physical 
modeling one or more of the following characteristics: fundamental-charge control 
physics for sheet charge in the active channel as controlled by the gate terminal voltage; 
average centroid position of the sheet charge within the active channel width; position 
of charge partitioning boundaries as a function fo gate, drain and source terminal 
voltages; bias dependence of linear channel mobility in the surface depleted regions; 
bias dependent of the velocity saturating electric field; saturated electron velocity; 
electrical conductance of the linear conductance region of the channel under the gate; 
electrical conductance within the source and drain access regions. 

14. The method as recited in claim 13, wherein step (c) includes the (e) 
adjusting the empirical terms of the semi-physical model to fit I-V characteristics 
against measured values. 
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15. The method as recited in claim 14, wherein the empirical terms of the semi- 
physical model are interactively adjusted to simultaneously fit the measured I- V and C- 
V characteristics. 

16. The method as recited in claim 15, wherein the empirical terms are fired for 
future use. 
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